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We discuss the role of higher order couphngs in conjunction with kaon condensation using recent 
versions of relativistic mean field models. We focus on an interaction (G2) in which all parameters 
are obtained by fitting finite nuclear data and successfully applied to reproduce variety of nuclear 
properties. Our results show that the higher order couplings play a significant role at higher den- 
sities where kaons dominate the behavior of equation of state. We compare our results with other 
interactions (NLl, NL3, Gl and FSUGold) and show that the new couplings bring down the mass of 
neutron star (NS), which is further reduced in the presence of kaons to yield results consistent with 
present observational constraints. We show that the composition of NS vary with the parameter 
sets. 



PACS numbers: 26.60.-c, 26.60.Kp, 13.75.Jz, 97.60.Jd 



I. INTRODUCTION 

Neutron stars (NS) provide us with opportunities to 
probe the properties of matter at extremely high densi- 
ties, and have proven to be fantastic test bodies for the- 
ories involving general relativity. In a broader perspec- 
tive NS provide access to the phase diagram of matter at 
extreme densities and temperatures, which is the basic 
for understanding very early Universe and several other 
astro-physical phenomena. The observational quantities 
of primary astrophysical interest are the maximum mass 
and the typical radius of a NS. Neuron stars are detected 
as pulsar -|-NS or pulsar -|-white-dwarf or X-ray binaries. 
Recent observations of pulsars and X-ray binaries suggest 
that the maximum mass of a NS lies between 1 — 2 Af© 
[Il-il, where Mq is the solar mass. 

To understand the observables of NS, various theoreti- 
cal models have been developed and can be grouped into 
three broad categories : nonrelativistic potential mod- 
els relativistic field theoretical models [ll[ and 
Dirac-Brueckner-Hartree-Fock models 12 . [Tst . In each 
of these models, the addition of hyperon 14| or kaon or 
pion [isl - fTsj or quarks [l6j or their combinations, will 
soften the equation of state (EoS) and hence lower the 
maximum mass of NS. In all of these models, coupling 
constants and unknown meson masses are treated as ef- 
fective parameters adjusted to fit the empirical quantities 
(saturation density po, binding energy E/A, compression 
modulus Koo , effective nucleon mass m* and asymmetry 
energy J) at nuclear saturation. The compression mod- 
ulus (Koo) defines the curvature of the EoS at saturation 
density and its value will be reflected in the high density 
behaviour (stiffness or softness) of the EoS. Thus if 00 will 
have a direct bearing on the maximum mass. Interest- 
ingly, the study of isoscalar giant resonances can reveal 
rich information about nuclear compressibility [l^ [20| . 
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The EoS and hence the NS radius can directly be linked 
to the neutron skin thickness in heavy nuclei (2l| - [23| . 
These strong links between the NS structure and finite 
nuclear properties prove that the knowledge in these two 
areas complement each other. 

One among the very-well tested model in finite nuclei 
is the nonlinear a — lo — p model, widely mentioned as rel- 
ativistic mean field (RMF) model [H. The RMF model 
constructed with idea of renormalizability, produces ac- 
curate results for spherical and deformed nuclei but it 
gives a very stiff EoS for infinite nuclear matter. Re- 
cently, inspired by effective field theory (EFT) Furnstahl, 
Serot and Tang [25| abandoned the idea of renormaliz- 
ability and extended the RMF theory, with the system- 
atic inclusion of new interactions by adding new terms to 
the model Lagrangian. This effective field theory moti- 
vated RMF (E-RMF) model calculations explained finite 
nuclei and nuclear matter with exactly same parameters 
and with significant accuracy in both the cases [2^ . This 
approach can be considered as a salient step towards a 
unified and accurate theory for finite nuclei as well as 
for infinite nuclear matter. Extension of this model with 
the inclusion of the kaon {K~) is the central interest of 
present work. 

In NS, if "-condensation is one of the several possi- 
ble transitions that could exist at high density. As the 
density of NS increases, the chemical potential of the 
negative electric charge {pe) also increases with the same 
rate of proton number density. Simultaneously, the ef- 
fective mass of an in-medium will decrease, due to 
the attractive interaction between K~ and nuclear mat- 
ter. Therefore at a particular density (when p^ is greater 
than energy of A'"), K~ replaces electron as neutralizing 
agent in charge-neutral matter. After the demonstration 
of interaction of K~ with the nuclear medium by Kaplan 
and Nelson [23|, this topic acquired enormous interest. 
Glendenning and Schaffner-Bielich [13, [3 : explained the 
interaction of A'~ by coupling them to meson fields using 
a minimal coupling. This approach was followed by sev- 
eral other theoretical groups using RMF models [28l - [30| 
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and has been adopted in this work as weh. 

In the present work we investigate the condensation 
from non-kaonic to kaonic phase, with different La- 
grangians. In the presence of such a transition, we study 
the effect of higher order couphngs on the EoS and hence 
on the NS properties. This work differs from many pre- 
vious works in the sense that aU the parameters used in 
E-RMF model were fitted to reproduce the observables 
of finite nuclei. Also, the E-RMF model is well tested 
throughout the nuclear chart by explaining several nu- 
clear properties [3l| - l36j . We show in this paper that this 
model explains the recent observations of NS as well. 

In the section following this introduction, we describe 
the Lagrangian and field equations, both for non-kaonic 
and kaonic phases. This is followed by the expressions 
for energy density and pressure, which define the EoS. 
In section Hill we discuss the parameters used in our cal- 
culation. Our results and discussions are presented in 
section HVl which is followed by the summary along with 
the conclusions drawn from present work. 



II. EFFECTIVE FIELD THEORY MOTIVATED 
RELATIVISTIC MEAN FIELD MODEL WITH 
KAONS 

In this section we briefly sketch the E-RMF model by 
presenting the model Lagrangian [2^ [2^ . We then show 
how the physical quantities, that will determine the com- 
position of the NS, can be obtained self consistently. 

The effective Lagrangian, obtained by curtailing terms 
irrelevant to nuclear matter in the E-RMF Lagrangian, 
can be written as 



.9crO- 



9 1 



2 2 

-mid 



1 ^ K^g^a 

2 3!m„ 



2 2 



(1) 



where the scalar, vector and isovector meson fields and 
the nucleon field are denoted by cr, V^, R^^ and respec- 
tively. TTio-, m^j, and nip are the corresponding meson 
masses and ?ti„(— nip) is the nucleon mass. The symbols 
ga, guj, gp, ^3, ^4, 771, 772, Vp and Co denote the various 
coupling constants. More details of the Lagrangian are 
explained explicitly in Ref. [2^ . 

Now we extend the E-RMF model by including a kaon 
(if ~)-nucleon interaction term. For simplicity, we choose 
that the kaon is coupled to the meson field with minimal 
coupling [T^]. In this way, interactions of mesons with 
nucleons and kaons are treated in the same footing. 

The Lagrangian for the kaon part reads 
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D*K*Di'K 
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where the vector fields arc coupled to kaons via the rela- 
tion 



D^^=df_,+ igujRVp, + igpKTj, ■ i?^. 



(3) 



and m*^ is the effective mass of kaon. The scalar field 
is coupled to kaons in a way analogous to the minimal 
coupling scheme of the vector fields: 



= niK ~ gaK<y, 



(4) 



where wik = 495 MeV. Note that in the mean field 
approximation, only the time components of the vector 
fields Vo contribute and charge conservation implies that 
only the third component in isospin-space of the isovector 
meson field i?o does not vanish. 

The dispersion relation for s-wave condensation {k — 
0), for K- is 



^A' = mx - gaKO- - gojxVo - gpuRo, 



(5) 



with K~ mesons having isospin projection —1/2. lok rep- 
resents the K~ energy and is linear in the meson field. 
One can note that with the increase in density, ujk de- 
creases. 

In the presence of K~ , equations of motion for the 
meson fields are 
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Cogly^ - guKPK, 



m^pRo = T^gpiPp - Pn) - ilp'-^m^RQ - gpKPK, (6) 



where ps is the scalar density given by 



7r)3 ^ 



{2iTf ^^J, -'^{k^+mfy/ 

and 7 is the spin-isospin degeneracy factor and is equal 
to 2 (for spin up and spin down). The nucleon effective 
mass is defined as in the standard Walecka model 

777* = 777„ - gacr, (8) 

and the K" density is given by 

PK = 2m},K*K = 2{l,k + guxVo + gpKRo)K*K. (9) 

The proton and neutron chemical potentials are 

/"p = \J ^/p + +guiVo + igpi?o, and 

l^n = \J kj,^ + 777*2 ^ g^Y^ _ ^gpRo- (10) 
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For a NS, in the absence of neutrino trapping, the con- 
servation of baryon and electron chemical potentials leads 

to m 



(11) 



and 



q = Pp - Pe - Pp. - Pk, (12) 

where the first two constraints ensure the chemical equi- 
librium and the last one specifies the total charge which 
vanishes while imposing the charge neutrality. The total 
baryon density is 

P^Pp + Pn (13) 

The energy density can be written as 

e = eN + eK, (14) 
where e^v is the energy density due to nucleons given by 



B. Kaonic phase (n, p, e , ^ , K ) 

With the solution of non-kaonic phase in hand, from 
Eq. ([5]) we can calculate kaon energy which keeps decreas- 
ing as we increase density, while fi^. increases. When the 
condition loh = fi^ is first achieved, the kaon will oc- 
cupy a small fraction of the total volume and the charge 
density corresponding to kaonic phase, q = 0. We can 
calculate a, Vb, Ro, kfp, kfn, kfe, kf^ and px by using 
Eqs. dH), ((ni) dm) and with the condition ujk = Pe, 
for any chosen baryon density. After getting these solu- 
tions, we can calculate energy density and pressure for 
kaonic phase, using Eqs. and ([T71) . 
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and the energy density contributed by kaons is 

ex = 2m*^K*K ^ m*jiPK. (16) 

Unlike the energy density, pressure is not directly affected 
by the inclusion of K~ , but the inclusion of affects 
the fields and hence the pressure which reads 

1 
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Here / stands for the leptons (e^,/i^). 



A. Non-kaonic phase (n, p, e , p ) 



(17) 



III. 



CHOICE OF PARAMETERS 



In the present work we have chosen five sets of pa- 
rameters namely NLl [H, NL3 [H, Gl, G2 [13 and 
FSUGold [12 . The first two parameter sets correspond 
to the standard RMF model with very different compress- 
ibilities [K^{NL1) = 212 MeV and A'oo(NL3) = 271.76 
MeV]. In the case of FSUGold, there are two more cou- 
pling parameters A and Co in comparison to NLl and 
NL3 parameter sets. These parameters represent the 
strength of the self-interaction of vector field (Co) and 
the isoscalar -isovector mixing (A). The interactions Gl 
and G2 correspond to the complete E-RMF Lagrangian 
discussed in the earlier section. One can note that in 
comparison to NLl and NL3, Gl and G2 have four more 
coupling constants (771, 772, ?7p, and Co)- AH the coupling 
constants are obtained by fit to several properties of fi- 
nite nuclei It is worth to mention that in this case 
no parameter is treated as adjustable and the fitting does 
not involve any observation at densities above the satu- 
ration value. Since the expectation value of the Rq field 
is typically an order of magnitude smaller than that of 
the Vq field, in the E-RMF model, the nonlinear Rq cou- 
plings were retained only through third order j25i] . How- 
ever, it has been shown that the isoscalar-isovector mix- 
ing (A(g^Vb)^(gpi?o)^) is useful to modify the neutron 
radius in heavy nuclei while making very small changes 
to the proton radius and the binding energy [2l|, . 

In the effective Lagrangian approach adopted here, 
knowledge of two distinct sets of coupling constants — 
one parametrizing the nucleon-nuclcon interaction and 
one parametrizing the kaon-nucleon interactions — is re- 
quired for numerical computations. We discuss each of 
these in turn. 



In the non-kaonic phase, vanishing charge density im- 
plies q = with Pk = 0. We can calculate a, Vb, i?o, fc/p, 
kfn, kfe, and fc/^ by using Eqs. ©, JlTl), (US) and ([T^. 
at the chosen baryon density. When we get the converged 
solution for the above-listed quantities, the energy den- 
sity and pressure can be computed from Eqs. (jl4p and 



A. Nucleon coupling constants 

The symbols g^, gp, K3, K4, 771, 772, 77^, Co and A 
denote the nucleon coupling constants. At times these 
constants bear different values in literature for the same 
parameter set. For example, for the set NLl different 
numbers are quoted in Refs. (38j and [39| but they yield 
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TABLE I: Parameters and the saturation properties for NLljsll, NLSfsll, Gl, G2[3l, and FSUGold[2|^ The parameters g^, 
g^, gp, K3 and K4 are calculated from the given saturation properties using relations suggested in Ref. [33] with the exception 
of FSUGold where all the coupling constants are taken from Ref. [2^ . 





NLl 


NL3 


Gl 


G2 


FSUGold 


m„ (MeV) 


938 


939 


939 


939 


939 


(MeV) 


492.25 


508.194 


507.06 


520.206 


491.5 


rriu, (MeV) 


795.36 


782.501 


782 


782 


783 


rup (MeV) 


763 


763 


770 


770 


763 


9t 


10.0730 


10.1756 


9.8749 


10.5088 


10.5924 




13.1917 


12.7885 


12.1270 


12.7864 


14.3020 


gp 


9.8553 


8.9849 


8.7886 


9.5108 


11.7673 




1.8324 


1.4841 


2.2075 


3.2376 


0.6194 




-7.7099 


-5.6596 


-10.0808 


0.6939 


9.7466 


Vi 








0.071 


0.65 





rj2 








-0.962 


0.11 





Vp 


n 





— fl 979 




n 


Co 








3.5249 


2.642 


0.06 


A 














0.03 


Po (fm-^*) 


0.154 


0.148 


0.153 


0.153 


0.148 


(E/A)(MeV) 


-16.43 


-16.299 


-16.14 


-16.07 


-16.3 


A'oo(MeV) 


212 


271.76 


215 


215 


230 


J(MeV) 


43.6 


37.4 


38.5 


36.4 


32.59 


m*„/mn 


0.571 


0.6 


0.634 


0.664 


0.609" 



"m* for FSUGold was calculated from the coupling constants. 



same saturation properties. Among the above-mentioned 
coupling constants, Qa, g^, gp, ^3 and K4 can be written 
algebraically in terms of einpirical quantities: po, E/A, 
Koo, J, 'Ti* and vice versa [37[. Using these relations wc 
have calculated g^, g^, gp, and K4 in the case of NLl 
m, NL3 Hi], Gl and G2 Hi. These calculated param- 
eters, and other parameters taken from corresponding 
references are listed in Table HI In case of FSUGold, all 
the coupling constants are taken from Ref. [2^ . 



B. Kaon coupling constants 

In order to investigate the effect of kaons on the high 
density matter, the kaon-nucleon coupling constant has 
to be specified. The laboratory experiments give informa- 
tion only about kaon-nucleon interaction in free space. In 
this work we mainly focus on the NS (densities ^ po), 
and therefore the kaon-nucleon interaction determined 
from experiment need not be appropriate for our calcu- 
lations. The interaction of omega and rho mesons with 
kaon {g^jK and gpK) can be determined using a simple 
quark and isospin counting argument given by, 

guK = ^gu: and gpK = ^gp ■ (18) 

We can specify the interaction of sigma meson with kaon 
{g<yK) using its relation with the optical potential of a 
single kaon in infinite matter (Uk)' 

Uk = -gaK<y{po) - guiKVoipo), (19) 



where typically we have -80 MeV < Uk < -180 MeV 

Em. 



IV. RESULTS AND DISCUSSIONS 

The quality of different interactions in explaining sym- 
metric matter and pure neutron matter properties at 
higher densities than the normal nuclear densities has 
been discussed in Refs. [H, |42|, by comparing the cal- 
culated pressure with the experimental data [43| . This 
comparison clearly suggests a softer EoS for symmetric 
matter, which could be obtained only with higher order 
couplings in the Lagrangian. The pure neutron matter 
data, though model dependent, also was found to favour 
a softer EoS. One important observation highlighted in 
Ref. [2^, is that the softness in EoS from Gl, G2 and 
FSUGold is mainly due to the new couplings and not due 
to the difference in compressibility as usually perceived. 
For example, the compressibility from NLl (^^00=212) 
and Gl, G2 {Kao=215) are very similar but the EoS at 
higher densities arc completely different. NLl (_ftroo=212) 
and NL3 (A'oo=271.76) yield different compressibilities 
but their EoS are more similar. The role of higher or- 
der couplings in softening the EoS has been discussed in 
Ref. [il for the case of FSUGold. EoS from FSUGold is 
softer than that of G2 because of the large and positive 
K4 value as well as the introduction of isoscalar-isovector 
coupling (A). In fact it has been clearly demonstrated 
[21I . li^ l that A softens the symmetry energy consider- 
ably. 
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FIG. 1: (Color online) The density dependence of the 
energy (ojk) in NS matter for different optical potentials (JJk 
in MeV) calculated with G2 parameter set. The point at 
which electron chemical potential (/ie) intersects ujk defines 
the onset of Jf" -condensation. 



The success of Gl and G2 in explaining the high den- 
sity EoS constrained by phenomenological flow analysis 
[43| . has been carried forward in explaining NS prop- 
erties as well [2^. FSUGold also has been successful 
in explaining both these features [22|, |42j. Henceforth 
we discuss the role of new parameters of extended RMF 
models on the onset and effect of kaon condensation in 
NS. Among the different parameter sets considered here 
(Table m G2 and FSUGold have a positive quartic scalar 
self-coupling (K4), which is more meaningful than a neg- 
ative one [2a, |4^. For this reason we prefer the set G2, 
when we have to choose between Gl and G2. 



A. Results from G2 

In Fig. [11 we have presented the K~ energy (a;^) 
for different optical potentials, as a function of density 
along with the variation of electron chemical potential 
(/ie). The ojk decreases with increasing density [Eq. ([5])] 
while /ie increases. When u:k is lower than /te, K~ are 
favoured (due to attraction between K~ and nucleon) to 
replace electrons while contributing to the charge neu- 
trality. From Fig. [U it is evident that the onset of K~- 
condensation is strongly modified by the strength of the 
kaon optical potential (Uk)- Calculations based on chi- 
ral models [41| suggest that a value Uk ~ —120 MeV is 
more appropriate for NS. Hence in most of the further 
discussions we choose Uk = — 120 MeV and look at the 
other dependencies for the EoS of NS. 

The a, lu, and p fields calculated for NS at Uk — —120 
MeV are presented in Fig. [5] along with the /te. The 
density at which K~ starts to contribute can be read 
from the point where the /ie shows a sharp kink (^ 0.6 
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FIG. 2: (Color online) The density dependence of the scalar 
{gaKu), vector {qukVo), and iso-vector (gpKRo) fields in the 
NS matter inclusive of kaonic phase, calculated with G2 pa- 
rameter set. 
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FIG. 3: (Color online) The relative population of hadrons and 
leptons in NS as a function of baryon density calculated with 
G2 parameter set. The calculations done without considering 
kaons are represented by the small-dotted lines. 



fm^'^). Though the corresponding kinks in the fields are 
very subtle, any small change in the field results in a 
significant change in energy density and pressure. The 
presence of K" alters the proton-neutron ratio and hence 
the contribution from the p field is enhanced. This is due 
to the fact that the processes like 7i p + K~ are more 
energetically favoured than processes like n — > p + 
at higher densities where K~ energy is decreasing and 
/ie is increasing [See Fig. [Tj. From Fig. [5] we can also 
see that the K~ has very less infiuence on the a field. 
In other words, for the parameters considered here, the 
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FIG. 4: (Color online) Pressure in NS matter versus baryon 
density calculated using difTerent optical potentials and G2 
parameter set. 



attraction caused by is weak in comparison with that 
of a mesons. 

The appearance of K~ at higher densities and its role 
in population (number density) of different particles are 
presented in Fig. [31 As soon as iir~-condensation sets 
in, the population of K~ rapidly increases with density 
and hence allows the presence of more protons. Inter- 
estingly, at very high densities the number of protons 
exceed that of the neutrons, which is not the case if 
were not present. This can be seen from the correspond- 
ing deviation from the non-kaonic matter trends. Apart 
from the arguments given earlier, the reason for p-K~ 
pairs being preferred to neutrons is quantitatively very 
well explained by Glendenning and S chaff ner-Bielich [131 ■ 
However, the symmetry term in the energy obviously pre- 
fer a symmetric matter and hence hinders the protons to 
be populated much more than neutrons. All the above- 
mentioned effects have a significant role to play in deter- 
mining the EoS which is discussed in the following text. 

Our results for the pressure calculated with G2 param- 
eter set and with different optical potentials for K~ are 
presented in Fig. S) On first sight one can appreciate 
the strong influence of Uk in softening the EoS, which 
is similar to the results of many previous works. Also 
for lower values oi Uk, the graph suggests a second or- 
der phase transition from a non-kaonic to kaonic phase. 
Only with a very high value of UkQ^ —160 MeV) one can 
have a first order phase transition. This is rather consis- 
tent with the observation in Ref. [2^ where higher order 
interactions are included using the lowest order chiral La- 
grangian. However, the role of more general higher order 
operators was considered to be an open question which 
is answered in the present work. One may expect some 
changes in the onset of if "-condensation if we consider 
the presence of a mixed phase (of kaonic and non-kaonic 
phases) [l7| . The presence of such a mixed phase is nor- 
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FIG. 5: (Color online) The electron chemical potential versus 
the neutron chemical potential in NS matter calculated with 
different parameter sets. The kinks represent transition from 
non-kaonic to kaonic phase. 



mally possible in first order ^('"-condensation where one 
gets a dip in the pressure like the one for Uk ^ —160 
MeV as shown in Fig. |4l This argument implies that 
the presence of mixed phase is not favoured in our case 
for Uk < —160 MeV. The main observation from Fig. 2] 
is that even in the presence of higher order couplings, 
the presence of K~ dramatically softens the EoS and the 
softness is proportional to Uk- 



B. Comparison between different interactions 

So far we have discussed our results with the G2 pa- 
rameter set and in this section we analyze how our re- 
sults compare with different interactions considered in 
this work. In Fig. [5] we show the chemical potentials 
which indicate the onset of if "-condensation for differ- 
ent parameters. The maximum value of /Ze i^J'T°'^) 
notes the point at which the transition happens and is 
seen to be larger for Gl and G2. ^^"'^ is almost same for 
NLl and FSUGold and the least value is for NL3. /u™"^ 
also quantifies the number density of electrons which falls 
sharply when K~ starts to appear [See Fig. Both 
^max ^Yic values at which the transition happens 
(/i^), strongly depend on the interaction, /i^ for G2 and 
FSUGold are more closer resulting from a closer EoS. A 
sharper kink in Fig. [51 leading to a situation where there 
could be more than one value of fie for a given /x„, would 
have indicated a first order phase transition. It is clear 
that for Uk = —120 MeV, all the interactions lead only 
to a second order phase transition where a mixed phase 
cannot appear. 

In Fig.[S]we compare the number densities arising from 
the calculations with Gl and FSUGold, which can be 
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FIG. 6: (Color online) The relative population of hadrons and 
leptons in NS matter as function of baryon density calculated 
with (a) Gl and (b) FSUGold parameter sets. The calcula- 
tions done without considering kaons are represented by the 
small-dotted lines. 



compared with Fig. [3] as well. The results for Gl arc 
very similar to that of G2 except for the early onset of 
if^-condensation in Gl. Between Figs.|6](a) and (b), we 
can see a significant difference in the relative population 
between protons and neutrons at higher densities. This 
is mainly due to the interplay between uk and which 
are plotted in Fig. [71 The decrease in lok or increase in 
/ie favours more of K~ and hence protons, luk varies 
linearly with the meson fields as given in Eq. ([S]) and 
should reflect the stiffness of the fields and hence that of 
the EoS. We can see in Fig. [7] that the u:k for NL3 is 
quite different from others due to the fact that the EoS 
of NL3 is quite stifFer than others [l^ HH . This stiffness 
arises from the vector potential which grows almost as 
a straight line [s^l- The quartic vector self-interaction 
brings down the vector potential and makes the equa- 
tion of state and hence ojk softer in other cases. The 
other quantity determining the onset of A'~-condensation 
is ^e{= fJ-n ~ fJ-p) which has terms similar to that of 
the symmetry energy and is dominated by the Rq field 
[Eq. PH)) ] . Thus the density dependence of symmetry en- 
ergy directly affects fi^, and is crucial in determining the 
onset of A'~-condcnsation. NL3, Gl and G2 yield sim- 
ilar symmetry energy (not shown here) and the unique 
isoscalar-isovector mixing {A{g^Vo)'^{gpRo)'^) in FSUG- 
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FIG. 7: (Color online) Density dependence of the electron 
chemical potential (/ie) and K~ energy {uk) for different pa- 
rameters. The point at which /le intersects ljk defines the 
onset of A"^ -condensation. 
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FIG. 8: (Color online) The EoS for non-kaonic and kaonic 
NS matter obtained from various parameter sets. The cal- 
culations done without considering kaons are represented by 
the small-dotted lines whereas the calculations with kaons are 
represented by different line patterns as given in the legends. 
The solid circles correspond to the values at the center of 
maximum mass NS. Solid squares represent the observational 
extraction [46| . however not uniquely constrained [i^ . 

old suppresses the i?o field which leads to a softer sym- 
metry energy and hence lesser /.tg. 

If we compare between G2 and FSUGold, /ie(G2) 
> /ie (FSUGold) and with similar w/f, kaons are more 
favoured for larger /ie . Lesser kaons in FSUGold leads to 
lesser number of protons at higher densities in Fig.[S](b). 
In Fig. [7] one can see that the //e for NL3 and Gl are 
almost same till the transition point. However due to 
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FIG. 9: (Color online) The mass-radius relation for non- 
kaonic and kaonic phases using different parameter sets. The 
calculations done without considering kaons are represented 
by the small-dotted lines. The solid circles represent the max- 
imum mass in every case. Mass is given in units of solar mass 
Mq. Solid squares (t^k = R) and open triangles (r^u 3> -R) 
represent the observational constraints where r^h is the 
photospheric radius. The shaded region correspond to the 
recent observation of 1.97 ± O.O4M0 star 



the large difference in lok , which falls very sharply in the 
case of NL3, A'^-condensation occurs at a smaller den- 
sity in NL3 when compared with Gl. These features can 
be seen from the EoS presented in Fig. [5] as well. 

Figure [S] shows the pressure calculated with the differ- 
ent interactions, with and without the inclusion of . 
Similar to the case of symmetric and pure neutron matter 
(Hill, the higher order couplings in Gl, G2 and FSUG- 
old lead to a softer EoS. All the interactions suggest a 
major change in pressure when we include K~ . The well 
known feature of K~ making the EoS softer, can also 
be seen clearly in Fig. [51 For both kaonic and non-kaonic 
phases, the difference between pressure obtained from Gl 
and G2 changes with density. Around p ^ 0.5 fm~^, the 
difference is maximum and it decreases as the density in- 
crease. This is due to the interplay between the terms 
with higher order couplings {ks, H4,r]2) which can give 
negative contribution to pressure. Between Gl and G2, 
the sign of the coupling constants K4, 772, and rjp differ. 
Among these three couplings, the term comprising rjp is 
weaker. In the case of Gl, due to the change in sign of K4, 
the combined effect of K3 and K4 on pressure is negligible. 
Hence in Gl, 772 dominates at higher densities in reducing 
the pressure. In case of G2, the contribution from terms 
having K4 and r]2 is negligible and K3 dominates at higher 
densities in reducing the pressure. Thus in Gl and G2, 
772 (quartic scalar- vector cross-interaction) , and K3 (cubic 
scalar self-interaction) respectively, rule the suppression 
of pressure at very high densities. 772 being associated 
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FIG. 10: (Color online) The mass-radius relation for NS calcu- 
lated using different optical potentials and with G2 parameter 
set. Mass is in units of solar mass Mq and Uk are in MeV. 
The solid circles represent the maximum mass in every case. 



with the vector self-interaction is much stronger at high 
densities and hence will be more effective than K3 in re- 
ducing the pressure. Couplings defined by 7^2 and K3 are 
of different order and will naturally yield results of dif- 
ferent curvature which leads to the varying difference in 
pressure between Gl and G2. The difference in pressure 
from G2 and FSUGold is well reduced in the presence 
of K~, even at higher energies. Interestingly, G2 with 
K~ yields same pressure as that of FSUGold without 
and in a broad energy range these results cease to 
differ. However, the corresponding energy densities are 
very different (not shown here) and hence the resulting 
maximum mass NS have different central densities (rep- 
resented by solid circles in Fig. The EoS around the 
region of central density is more dominant in determin- 
ing the properties of NS. This region is quite different 
for different parameters. Also one can see that the cen- 
tral density increases with softer EoS, and the presence 
of K~ decreases the central density and pressure. 

The change in central density (Apc) due to K~ de- 
pends mostly on the following three quantities, viz.. 



(i) density at which K 
with Apc oc 1/pg, 



-condensation sets in [p^ ) 



(ii) the magnitude of central density with Apc oc pc 
and 
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FIG. 11: (Color online) The composition of maximum mass 
NS as a function of radial distance calculated using FSUGold 
and G2 parameter sets. 



(iii) the stiffness of EoS at pc, with stifFcr EoS leading 
to larger Apc- 

The interplay of the above quantities gives arise to vary- 
ing Apc for different parameters. For example, between 
Gl and G2, quantities (ii) and (iii) dominate (i), result- 
ing in a larger Apc in case of G2. These effeets should 
naturally be seen in the calculation of maximum mass 
and radius of NS. 

We can obtain mass-radius relation for NS by solving 
the well-known Tolman-Oppenheimer-VolkofF (TOV) 
equations [U 13 ■ The results for mass-radius relation 
in NS are given in Fig. |9] where the calculations are done 
using different interactions, without and with the K~ 
(Uk = -120 MeV). The RMF models NL3 and NLl 
suggest very large and massive NS. In these cases though 
the maximum mass is reduced considerably (from 2.8 
M0 to ^ 2.4 M0) in the presence of kaons, these are 
overshooting the observational constraints for NS [H-Q. 
Though Gl and G2 are from the same E-RMF model 
with same terms in the Lagrangian, their results for NS 
arc quite different with Gl suggesting a larger and heav- 
ier NS. As explained earlier, Gl and G2 yield EoS of 
different nature and have varying role of K~ . The result 
for G2 and FSUGold are the closest ones. Very interest- 
ingly, though their difference is larger when K~ are not 
considered, the presence of K~ brings these two results 



closer. This is the direct implication of the same feature 
we observed in the case pressure (Fig.|8|). The G2 inter- 
action allows more K~ to be present (due to early on- 
set of _fsr~-condensation) and hence makes the EoS more 
softer which eventually get closer to that of FSUGold. 
The observed reduction of maximum mass due to kaons, 
in case of NL3 and FSUGold is consistent with the val- 
ues quoted in Ref. [2^. Overall, in the presence of K~, 
the mass decreases and the radius corresponding to max- 
imum mass increases. All the observational constraints 
shown in Fig. IH] favor the models with higher order cou- 
plings but most of these constrains cover a broad range 
which may not be precise enough to ascertain the pres- 
ence of kaons. 

In Fig. [TOl we show the dependence of NS mass-radius 
on Uk with G2 parameter set. We can see clearly that 
the stronger Uk increases the contribution of K~ which 
results in softer EoS which in turn yields a smaller and 
lighter NS. The typical hook shape of the mass-radius 
curve is maintained till Uk ^ —160 MeV and beyond 
that the curve takes a straight dip which is a character- 
istic of first order i^^-condensation. It is interesting to 
note that without K~ , the maximum mass given by G2 
(I.95M0) is closer to the one suggested by Demorest et 
al. (1.97 ± O.O4M0) with a radius of 11.03 km. Gl 
with Uk = —140 MeV (graph not shown here) suggests 
a maximum mass of 1.97M0 with a radius of 12.45 km. 

Figure [Tl] shows the number density versus the radial 
distance from the center of NS, calculated with G2 and 
FSUGold parameter sets, when the NS is having the max- 
imum mass (1.68Af0 for FSUGold, and 1.82AfQ for G2 
with Uk = —120 MeV). The behaviour of number den- 
sities shown in Fig. [11] is directly reflecting the patterns 
shown for relative populations plotted against density in 
Figs. [3] and [6] (b). The difference between the radii of 
kaonic matter suggested by G2 and FSUGold is almost 
one km. Another striking difference in the composition 
of NS is the asymmetry of the core. G2 suggests a nearly 
symmetric matter at the core whereas FSUGold indicates 
highly neutron-rich core. We have seen that this trend 
remains same even at higher values ol Uk (not shown 
here). It will be interesting to see how the presence of 
hyperons affect the overall scenario [49l |. 



V. SUMMARY 

In the present work, we have studied the role of higher 
order couplings of extended RMF models on the onset 
and effect of kaon {K~) condensation in neutron stars 
(NS). We have calculated the NS properties with the suc- 
cessful RMF parameters NLl and NL3, the E-RMF pa- 
rameters Gl and G2 (four additional couplings to RMF), 
and the FSUGold parameters (two additional couplings 
to RMF). 

In extended RMF models, with most common values of 
kaon optical potential {Uk —160 MeV), the transition 
from a non-kaonic phase to kaonic phase in NS has a 
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character of second order which rules out the possibihty 
of having a mixed phase of kaonic and non-kaonic matter. 
This is consistent with the observation in Ref. [1^ where 
higher order interactions are included using the lowest 
order chiral Lagrangian. The role of more general higher 
order operators was considered to be an open question 
and wc have shown that they lead to a second order phase 
transition. This justifies the neglect of mixed phase in our 
calculations mostly done with Uk = — 120 MeV f4l| . 

Apart from the usual dependence on Uk reported else- 
where [l3, H^, the onset of kaon condensation in NS 
strongly depends on the parameters of Lagrangian es- 
pecially the higher order couplings. This is due to the 
strong variation in density dependence of the K~ energy 
(ujk) and electron chemical potential (^e) whose inter- 
play determine the onset of ^"-condensation. Density 
dependence of ujk is similar to that of EoS and /ie varies 
in a way similar to symmetry energy. So, any change 
in the density dependence of EoS or that of symmetry 
energy will affect the onset as well as the effect oi K^- 
condensation. Without higher order couplings, NLl and 
NL3 have stiffcr EoS than others which leads to a stiffcr 
ujk and hence a early onset of iC~-condensation. FSUG- 
old comprises an unique isoscalar-isovector mixing which 
leads to a softer symmetry energy and hence a delayed 
onset of _ftr~-condensation. 

The central density of NS increases with softer EoS, 



and the presence of K~ decreases the central density and 
pressure. The change in central density due to the pres- 
ence of K~ is different for different parameters due to the 
interplay between (i) the density defining onset of K~- 
condensation, (ii) the central density (pc) itself, and (iii) 
the stiffness of EoS at pc- Due to this, the impact of 
is more pronounced in G2 than in Gl and is weakest in 
FSUGold. All these effects arc strongly refiected in the 
calculation of mass-radius relation of NS. The NS sug- 
gested by models without higher order couplings (NL3 
and NLl) contradict the observational constraints [ll-Q 
even with the inclusion of K~. The higher order cou- 
plings play a dominant role (than kaons) in bringing the 
mass and radius of NS within observed limits. Different 
parameter sets lead to different concentration of kaons 
(which is known to affect the population of protons [T7| ) 
in NS, and hence lead to different asymmetries in the core 
of NS. We have shown that Gl and G2 suggest a symmet- 
ric core whereas FSUGold suggests a neutron rich core 
due to lesser amount of K~ caused by an unique coupling 
which softens the symmetry energy. 

We conclude that the extended RMF models, which 
are quite successful in explaining several properties of 
finite nuclei, suggest a strong influence of higher order 
couplings and kaons in NS whose properties are within 
the present observational constraints [lH6| . 
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